


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971-03 


An evaluation of a numerical water elevation 
and tidal current prediction model applied to 
Monterey Bay. 


Lazanoff, Sheldon Mark 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/15701 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


(8 DUDLEY research materials and institutional publications created by the NPS community. 
«ist Spe Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 


INN KNOX appointed — and published -- scholarly author. 

| LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





AN EVALUATION OF A NUMERICAL WATER ELEVATION 
AND TIDAL CURRENT PREDICTION MODEL 
APPLIED TO MONTEREY BAY 


by 


Sheldon Mark Lazanoff 





RES ER em ae RR a 


United States 
Naval Posigraduate Scncol 








cD Saye CURRENT PREDICTION MODEL 
APP GED TO OMONTERZY BAY 


DY 


oO 


4 
AN HEVALUATION OF A NUMERICAL WATER ELEVATION 
Sheldon Mark Lazanoff 
I 


bd ph, * 
oe 6 eee eee op 
t-. : 


a. Advisor: Eee nar inieoln 


March 1971 








Approved for public xckease; aistrioution wilimcted. 








An Evaluation of a Numerica_ Water Elevation 
and Tidal Current Prediction Model 
Applied to Monterey Bay 


by 


Sie laeieiarky Lananort 
Civilian, United States Naval Oceanographic Office 
Bowe Pennsyivania State University, 1963 


Submitted in partial fulfillment of the 
requirements for the degress of 


MASTER OF SCIENCE IN OCEANOGRAPHY 
from the 


NAVAL POSTGRADUATE SCHOOL 
March l9O7L 





LIBRARY 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY. CALIF, 93940 


ABSTRACT 


The Hansen Hydrodynamical - Numerical model was evaluated 
for Monterey Bay with actual field data. Tides and winds axe 
the principal driving forces of the Hansen model. Analysis 
of the field data indicated that the principal Grivang, £OLrCe 
of the circulation in the bay was the oceanic currents and 
not the tides and winds. The tidal heights and phases and 
current directions were calculated correctly by the model, 
but the calculated current speeds were an order of magnitude 
too large. The inaccuracy of the current speeds was attri- 
buted to the inaccurate calculations of the currents along 
the open boundary and the large bathymetric gradients of the 


Monterey Submarine Canyon. 
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I. INTRODUCTION 


A. REVIEW 

During the normal course of events, tidal circulation in 
coaStal areas, estuaries and eine affects a diverse number of 
paramenters including: marine life in the inter-tidal and 
near shore zones; the flow of man-made sewage and other pollu- 
tants dumped into the water; off-shore structures, channel 
mmecaging and other engineering projects; and, of course, ship 
traffic. Furthermore, in areas where intense storms occur, 
the wind can create storm surges which combined with the tides 
can kill and injure people and destroy property on the sur- 
rounding low lands. Before the advent of high speed digitai 
computers, the prediction of tidal heights and currents was 
difficult and cumbersome. The tidal circulation in deep water 
tmpbimarily due to the attraction by the moon and sun; however, 
the tides in shallow water are also strongly affected by the 
geography and bathymetry of the area. 

The techniques for tide medic edon developed by the US 
Coast and Geodetic Survey (USCGS) and the British Admiralty, 
among other international agencies, in the late nineteenth 
century have become standard throuchout the world. The pro- 
cedure requires tidal height and current data measured for 
at least 29 days which covers the principal lunar constituents, 
preferably for a year which covers the principal solar con- 


stituents and ideally for 18.6 years which covers all possible 


NZ, 





tidal constituents. The tidal constituents, or harmonic 
coefficients, are then used to predict the tidal elevation. 
The accuracy of the predictions increase with the length of 
the record. The data are analyzed by harmonic techniques and 
tidal constituents are produced. In the past the tidal con- 
stituents were used as input to large tide prediction machines 
which produced hourly heights and tidal currents (Schureman, 
1958). The machine needed at least twelve hours to produce 
yearly predictions for one station. One serious problem with 
this technique is that if tide data are not available for a 
location cf interest, then the predictions have to be inter- 
polated between known stations as shown in the USCGS Tide 
Tables. The interpolated tide values can be inaccurate and 
more Aeeameneal than helpful for real time operational 
planning. The U.S. Navy became acutely aware of this problem 
when using French tide tables for riverine and coastal oper- 
ations in Vietnam. 

The necessary equations of motion needed to describe the 
tides were developed by Bernoulli and Laplace in the 18th 
century. However, because the equations are non-linear, 
Simplfying assumptions had to be made and tedious analytical 
methods were used to obtain solutions. Analytical methods 
are not practical for real time predictions. By using 
computers these equations can be solved in less time by more 
sophisticated-methods. For example, yearly tide predictions 
for one station can now be calculated on a computer in 1.7 
minutes (Pore and Cummings, 1967) compared to twelve hours 


for the mechanical tide prediction machine. Several 
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computerized mathematical models using different numerical 


—_ eS 





schemes and boundary conditions have been developed to predict 


: 2 “ . _— 


ee 


assorted circulation parameters in a two dimensional field 


pate 
— 


rather than for one point. An operational hydrodynemical 


a — ER rn TEP at oa = 


model to compute water elevations-and—currents has been 
developed vy Professor Walter Hansen, Universitat Hamburg. 
Ree two driving forces in the model are the tides at the 
open boundary and the wind stress acting at the surface ove. 
the entire grid. The Hansen Model was initially begun in 
1938 (Laevastu and Stevens, 1969), culminating in a two 
dimensional model of the North Sea in 1952 (Mungall and 
Matthews, 1970). In 1966 the Nortk Sea model was programmed 
for the ccmputer by Jensen, Weywadt, and Jensen (1966). The 
predictions compared favorably with field observations. Some 
of the general characteristics of Hansen's two dimensional 
model are that density is uniform over depth; the water 
transport is averaged over depth; bottom and wind stresses 
are represented by quadratic formulas; Coriolis force is 
assumed to be constant over small geographical areas; the 
effect of atmospheric pressure is neglected; and there are 


no more than two open boundaries on the rectangular grid. 


poe ObeECLIVES 

Since June 1968 Fleet Numerical Weather Central (FNWC) 
and U.S. Naval Oceanographic Office (NAVOCEANO) personnel 
have modified and adapted the Jensen, et.al. program to other 
geographic areas such as the South China Sea, DaNang Bay, 


Gulf of Tonkin, Straits of Gibraltor, Gulf of Mexico and the 
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fo 


Chesapeake Bay. Even though Hansen himself has indicated 
(Hansen, 1966) that synoptic field data are extremely impor- 
tant in verifying the results of the computer calculations, 
the results for the above mentioned projects have only been 
compared to atlas et ane sparse field data, leaving some 
doubt as to the accuracy of the predictions. The U.S. Naval 
Oceanographic Office does oceanographic surveys in shallow 
water, semi-enclosed embayments throughout the world to fulfill 
specific requirements in supporting naval operations such as 
mine sweeping and diffusion studies. Because of ever changing 
world politics, oceanographic surveys cannot always be made 
in vital areas before they are neeced and then, circumstances 
Can prevent a proper survey from béing made in these areas. 
Such situations have occurred in Vung Tau and DaNang Bay in 
the Republic of Vietnam. Thus, the object of this report is 
to evaluate the accuracy of Hansen's mathematical model for a 
shallow water, semi-enclosed embayment; varying parameters 
such as boundary conditions and the horizontal eddy viscosity 
parameter to determine their effects on the predicted eleva- 
tions and currents and to determine the minimum field data 
requirements needed to satisfactorily verify tne computer 
results. 

‘For several reasons Monterey Bay appeared to be an excel- 
lent site to evaluate the model. Monterey Bay has similar 
geometrical features to many of the NAVOCEANO areas of 
interest. Unlike Southeast Asia areas there is no significant 
river run-off to influence the results (Lazanoff and Clarke, 


1970). Winds are fairly steady and predictable over a long 
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period of time allowing the model to be evaluated for two 
steady-staze conditions-summer caln. conditions and fall and 
winter storm conditions. Finally, Monterey Bay, which is 
bisected by a one of the deepest submarine canyons in the 
world allows the model to be testecl over a large change of 


depth in a short horizontal distance. 
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II. THE HANSEN HYDRODYNAMICAL MODEL 


A. MATHEMATICAL MODEL 


Although the Hansen two-dimensional model is detailed by 


moon, Ct. al., the following is @ review of the more impor- 


tant characteristics of the model. The coordinate system and 


surface and boundary conditions are shown in Figure 1. Some 


of the basic assumptions made in the analysis scheme are: 


I 
2 


3 


) 
) 


) 


the fluid is homogeneous and incompressible, 
the fluid is in hydrostatic equilibrium in the 
vertical direction, 

the geographical and the vertical variations 


of the Coriolis force is neglected. 


The equations used in the Hansen model are derived from 


the conservation of momentum and the conservation of mass 


which are respectively: the conservation of momentum 


where 


f 


ys 


t 


a 


Ou. du. e 

st a fa ttl) oe Sen. (i 
be My oe men paoee yogis (1) 
aes 


is the velocity component in the 2 direction, 

is time, 

is the Coriolis parameter, 

is the fluid density and, based on assumption (1), 
is constant over depth, 


is the pressure of the fluid, 


days 








h Inet hop, coy ts) 


a xe 


<-—_—__—_- Tt ————>} 
x 
Sei 


Figure 1. Hansen Coordinate System. 
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a 


a. represents the wind, bottom and lateral shear 
stresses combined, 

gs represents external body force which is here 
limited to the gravitational attraction of ‘the 
Samet llr, 


and the conservation of mass. 





00 Opus ; 
oS ++ aX, = qal4 27s (2) 


The density of sea water is Easuned constant by assumption 
(1). This assumption is a limitation of the model and was 
made because of computer memory and time limitations at the 
time of the development of the original model. A multi-layer 
model is being developed at Fleet Numerical Weather Central 
under the guidance of Professor Hansen and Dr. Taivo Laevastu. 
The model was not available for use in Monterey Bay. 

From assumption (2), equation (2) in the vertical direc- 


tion reduces to 


n 
p= | p(X.) g dx, + P. 


—h 
(3) 
= og (nth) + Py 
where Es is the atmospheric pressure. 
The boundary condition at the bottom, X3 = 0, is 
U.'n = 0 (4) 


ak 


where n is the unit normal to the bottom. 
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The boundary condition at the surface, X4 =nth is 


on on ae Ce eee 
ane oxa 0 eee >) 


ee 3 
where n is the water surface elevation and U3 is the velocity 
in the vertical direction. 

The baratropic atmospheric pressure term has been elimin- 
ated from the Monterey Bay model. Admittedly, the water 
height wili increase by one centimeter as the atmospheric 
pressure decreases or increases one millibar (Lazanoff, 1969); 
however, this fluctuation is relatively insignificant compared 
Peete irlUctEUations Caused by the wind stress (Jensen, et. al.) 
and over the distance of Monterey Bay the atmospheric pressure 
varies only slightly except during very stormy periods. 


; du. 
The convective terms, u. 1, has also been neglected 


dX, 

from the Hansen model. If a rather large value of 100 cm/s2c 
(about 2 knots) is assigned to both Ue and the change of u; 
with respect to — and the distance between grid points is 
900 meters (actual mesh length for Monterey Bay), then the 
convective term has a value in the order of Tomes SiS 
Under the same conditions, the local Becelenati on term will 
be of order one; the Coriolus term will be on the order of 
one cm/sec and the gravitational term will be at least on 
the order of fae ay eccer Thus, the convective term will 
at least two order of magnitude less than any of the other 
terms in the equation. 

By integrating the remainder of terms in equations fy) 


and (2) over depth, combining terms and defining the mean 
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velocity as 


ae eel ca 
rE | u, dX, =) ay2 (6) 
~h 
where H = nth and is the total water depth, the following 


equations arise: 





Z t “+ 
ou. Ce nt cae ne 
1 me (a) = ae ua do on 
Pt 1) fe U, v re: = oH Fe (7) 
ae 
J 
yea lire 
and 
d (HU. ) 
on ee. 
aE + aK, = 0 (8) 


where v is horizontal eddy diffusivity coefficient (lateral 


shear stress coefficient), 


wy 1s the wind stress in the x direction in dynes /cm 


tC 1s the bottom shear stress in the Xs direction im 
dynes /cm* 

Equations (7) and (8) are the equations used in the Hansen 
model. 


The wind stress term Cie ), aS used in Hansen's model, 
7 


is assumed to bea quadratic expression in windspeed. This 


expression is 


Zpeee a (9) 


Ww _ 2 
Ls or W, (Ww + W.) 


Zu 





a 


where 4 is the wind drag coefficient (dimensionless) and Ww. 
and W. are the wind components. 
The bottom stress term Ge ) 1S represented similiarly 


1 
to the wind stress term and is 


a 2 272 
Ty =pr U; (U + a (10) 


where r is the bottom friction coefficient (dimensionless) 
and 1s assumed constant. 

Both the bottom stress and wind’ stress terms are nonlinear 
and have been derived by empirical means. This formulation 
is an approximation and its applicébility to deep water is 
questionable as it was originally formulated for shallow 


water application. 


Eee COMPUTER MODEL 

A central differencing scheme iis used to program equations 
(7) and (8). The grid system is shown in Figure 2 (Laevastu 
and Stevens, 1970). A "leap frog" method is used whereby the 
grid is staggered in time and space (Mungall and Matthews, 
1970). N and M correspond respectively to the Xs and Xo 
coordinates. The Z points are the symbolic depths indicating 
land, sea and boundary points. 


The finite approximations (Ortiz, 1964) are: 


Geet) —(t) 2 E) Ce) . (et) (les) 
n,m) ~ Tovm) ~ 77” (Fu, ae 
(N+1,M) (N-1,M) 
Uy ce! se 2 Hen ee ; 
(N-1,M) ge) tN M-1) 2 eg Mai) | eM) 


Ze 








. 1 2 3 4 
Z AU 
IZ b 
2 
N ) 
3 
ff 
a 
x Z-point, e ZU~point, Oo ZV-point. 


Pigure 2. Model Grid System. 
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where t is the time step between two Z points and 2 is the 


mesh length between two Z points 


xX (t) 


Get2t) _ aft ) (t) L 
U = (1-9 ) + 2tT£Us - -g (12) 
(NM) i M) ey M) (N,M) Ne 
(t+T) n (t+) x W(t) 
f (N+1,M)_ “Dao oe 
* (NM) 
Similarly, 
Sal / 
(t) cia iG 
ae Seg cee - ~g 
2 (Nn, M) > 1 L C3) 
(ceca Zz n ttt) \ W(t) 
f (N,M+1) ~ "(N,M- 1 et ie 
(N ,M) 


The terms used in the previous equation can be defined as: 


Ce) 

H = hu (t-T) (eat) 
ee Ponmy * 1/4 h (1+1,M) ~ "(N-1 : an 
(t) = Aes T ) CcSt) 

at = hu, F172 (No Ma1y) = (N, M41) MEP 
2 (NM) (N ,M) 


on = 2rT ce } + o (t) }? oe? (lise 
al (N ,M) (N ,M) (N ,M) 1 (N,M) 


One = an (of | oe ra 2 eager (17) 
(N ,M) 2 
2 (N/M) saa (NM) 
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UF and Us are the averaged velocities for the four Zu and 


Zv points about a Z point. 


Ima j)(t) 


ott) = aus”? + 7 I + as (18) 
. (N ,M) (N ,M) : (N+1 ,M) (N-1 ,M) 
+ ihe? cue 
(N,M-1) (N,M+1) 
ate ce ees + ust) (19) 
U = aU + —— |U 2 
2 (NM) 2 (NM) 4 \ 2 etm) Nee) 
ust) 
(N ,M-1) (N ,M+1) 
Ge?) = (ey) ac (t) ee), 
Di) ~ my) ~~ “4 | enti,my) * "-1,m) (20) 
mete) (t) 


(N,M-1) * "(N,M+1) 


Alpha (a) is a function of the horizontal eddy diffusivity 
ceriremesi Length and time (0 < a < 1). This term will be 
discussed in more detail in a later section. 

If any of the above calculations occur at a boundary, 
then the values of n, Us and U, are Pacer at the actual point 
rather than from the surrounding points. 

The primary forcing function for this model is the tide 
at the open boundary points. When the model was applied to 


other areas, predicted tides were generated for each time 


step using the following equation (Schureman) : 


= + 
a St 


I ma 


tL AL cos fa, t + Vee) _ Kn! (21) 


n=1 
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where, Nn = height of the tide at any time 


A, = mean height of water level above datum line, 
AL = amplitude of the ae amplitude 
fn = factor for reducing amplitude A to year of 
predictions, 
= ey 5 , 
AL = speed of the n constituent, 
c = time calculated from some initial epoch such 
as beginning of year of predictions, 
TU.) = value of equilibrium of the ate constituent 
when t = 0O, 
: th . 
> epoch of the n constituent. 


Equation (21) is applied only at the open boundary. Since 
eeemoegecel ve sof this particular project is to evaluate 
Hansen's hydrodynamical model and not the prediction of 
astronomical tides, a different method was used to calculatsa 
the input at the open boundary points. A Fourier analysis 
was made of the actual tide records which are used in this 
study. The Fourier constituents were used in the following 
equation:- 


N 
ys 


= l]/2 A + BS Ot + Be SiN Ot (22) 
n Vy S oe : Oe, 


A=1 


whexe, Aor me and ae are the Fourier constituents, determined 
from the actual records 
N is the total number of constituents used, 


Jj is any given constituent, 


Oo = on and T is the length of time that spans the 


entire record. 
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The principal parts of the Hansen hydrodynamical model 
have now been detailed and the following section describes 


the application of the model to Monterey Bay. 


Gee MIONTEREY BAY 

Monterey Bay (Figure 3) is a semi-enclosed (one open bound- 
ary), approximately symmetrical embayment. Because Monterey 
Canyon very nearly bisects the bay, there is an unusually large 
variation of depth (400 fathoms) in a relatively small area 
(approximately 175.5 square miles). The depth of the canyon 
thalweg from the mouth of the bay to the shoreline is shown in 
Figure 4. The thalweg from the shoreline to six miles off <zhe 
coast has an extremely steep slope of 1:18. Two transverse 
profiles cf the bay are shown in Figure 5. The gradients of 
the canyon wall is steeper on the south side than on the north 
Side but both canyon walls are extremely steep. Lynch (197)) 
found that the canyon tends to divide the bay into two dis- 
tinct basins in which long period oscillations such as seicne 
SecUureimaependent Of each other. The results of the Hansen 


ee 


model appear to indicate that this also occurs for tides which 


eee terete tna AP AIP trtati, due 
ete nmanttete nie 


qe dumin 


peeeecnde. periods than seiches: 

The total circulation of Monterey Bay is affected by 
density gradients as well as tides, wind stress and air _ 
pressure gradients. All these mechanisms can be important 
and it is not really correct to separate the density gradients 
from the other parameters. Unfortunately, the complex com- 
puter models needed to combine these parameters would consume 


too much computer storage and running time to be practical 


and economically feasible on the present day computers. 


2) 





' $ 
wee 504 
PA oa Nee 
tanra CRUZ i 
OXY >> 7 Sor S 
e 6 
Sees B 
*~, 
~.. 
Ne, 
Sx. 
~ 


a Station ] 
ae yx Current Meter 
fC 40 Feet . 
f--, 12 August - 20 August 1970 | 
Station 3 =. 
Current Meter x On . > 
90 Feet a oy 
7 Nov. - 12 Nov. ay 



















4 
. 
~ ‘. 
, \ Su 
rs | 
/ 4 ™“ 
€ a ~S 
pe a 3 
a ee \ 
3 ¢ J c ' = 
ony ee : = as \ ‘> a 
ae a ee ‘ eS . = eee oe F 
2 ee .. 3 { Se vcs 7 aon ¢ im ee 
a =, bf al 1 ae Saco St . a ‘ 
ne me ‘1-7 = oa “ 
SK / a ~~ es - J & 
«7% poner Te i” aoe of ~ o 
~") ~ eS % ui 
~~... fgg: / / — =.= =~ N \ Ys 
\ / : / ; : 
a“ / Pig : 
.- | Ci gg ae as / 
“-) : NN ae 
. \ ne f 
: s ‘ { 
x i. ; 
‘ \ { , 
\ ' \ y & 
ee / ‘ : ‘ 7 
. r / : : 
i) ; i | 
re ‘ : ‘ 
: \., a ea ‘ | 
.° 6 oats. } pe : 
¢ nae — a ae 4 a } D 
> mee : ! 
} Za ; 
. i am y “ 
: : 2 : / 
ne ; Station 2 ! 
: . Z ’ / 
: \ , 
; ss | : Current Meter, 3e°<ot 
fs o : Le 90 & eeu 
we ‘ 7 Nov. - 12 Novy 
x" / aN / 
e 5 . é 
a : if 
i 
a 
: Copy of C.&G. S. 
® 
, Q Map 5402 
° ee 
e ¢ 
i if ooqmecan +s =m 20 FATHOMS 
; 50 rT 
u socoommmm ef OO 
ah ee0e0 ap eee? = i000 
j 
Pace ax 


Figure 3. Tide, Wind and Current Meter Stations. 
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mewever, tO complete the picture of the circulation of 
Monterey Bay, some mention should Fe made of the semi-perma- 
nent density currents that occur in the ocean adjacent to 
the bay. | 

Depending on the climatological season, the circulation 
in Monterey Bay appears to be strorgly influenced by the 
California or Oceanic current, Davidson current and upwelling. 
The California current can be as much as 420 miles in breadth 
and have an average speed of 0.5 knots (Carter and Kazmierczak, 
1968). The current flows south during September, October and 
a portion of November. During the months of November througn 
February, the Davidson current, whi.ch is approximately 40 
miles wide and has a maximum velocity of 0.44 knots, flows 
north. The upwelling period extends from March througn July. 
mieeer, et.al. claim that the upweiling current flows south 
in the open ocean and north along the coast. The upwelling 
current velocities average about 0.5 knots over depth with a 
Maximum speed of 0.6 knots at the surface. The influence of 
the Monterey Canyon on the oceanic currents is not really 
known and the currents created inside the bay do not neces- 
sarily flow in the same ese rien as the oceanic current. 

1. Grid 

The choosen grid size is 19 (ordinate) by 46 {abscissae) 

with a mesh length of 0.5 nautical miles (Figure 6). The 
bottom depths were obtained from C & GS Chart No. S40 Seco care 
1:50000) which was updated in June 1970. The grid was 
selected so that the open boundary coincides with the nautral 


opening of the bay and Monterey Canyon 1s adequately defined. 
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Figure 6. Hansen Grid for Monterey Bay. 
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2. Time Step 


The time step between grid points was determined by 


using the Courant-Friedrich-Lewy criterion which is 


De 72 
DE i/légh, ce 
where, DT = 1/2 time step, 
DL = 1/2 mesh length (1500) ft. 
eae Maximum bottom depth of the grid area (400 
fathoms ) 


In this case the maximum DT is 3 seconds. 

If a time step greater than the maximum allowed by the 
Courant-Friedrich-Lewy criteria is selected, then the results 
of the model rapidly diverge. The author has experimentally 
verified this when applying Hansen's model to the South China 
Sea and the Gulf of Mexico. 

It is interesting to note that when equation (23) is 
rearranged, then the formula is different only by a factor 


of /2 from the deep water wave celerity formula: 


DL « 
pr ** "9%max ei 
3. Tides 


Since the tides are the primary forcing function 
across the open boundary, they need to be accurately defined. 
Ideally the tides should be known at every point on the open 
boundary; however, this is impractical. The optimum situation 
for Monterey Bay would have been to measure the tides at 


three points along the open boundary - at both ends of the 
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boundary and in the canyon which bisects the open boundary. 
Unfortunately, the tides were only measured near the open 
boundary at Monterey Wharf No. 2 (Point A on Figure 3) which 
is 3.3 miles from the closest open doundary point at the 
Southern end of the grid and the Santa Cruz Municipal Pier 
(Point B on Figure 3) which 0.3 miles from the closest 
northern open boundary point. It had been planned to place 
a submerged tide guage in the Canyon but the necessary 
equipment could not be assembled before the field survey 
eommenced. 

Hourly tidal height records from Monterey and Santa 
Cruz are shown in figures (A-1) through (A-6) and for Moss 
Landing in figures (A-7) through (4-10). The Monterey 
records have been corrected to a datum level. The tide 
records fcr the other stations contain the correct phases and 
relative emplitudes. The relative amplitudes for the three 
eratelous are compared to each other in Table A-1l. The maximum 
difference between amplitudes of any-two-stations is less tran 
0.75 feet and for most cases there was no difference at all. 
The Santa Cruz and Moss Landing records compared more favor- 
ably with each other than they did with the Monterey Wharf 
records. The fact that the tides at Santa Cruz and Moss 
Landing were measured by the same type of instrument while 
the tides at Monterey were measured by a different type of 
recorder may explain part of the difference. Since the tidal 
phases and amplitudes at Santa Cruz and Monterey were very 
nearly equal and there was no knowledge of the tidal behavior 


in deep water, only the Santa Cruz data were Fourier analyzed. 
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The Santa Cruz constituents were used at every open boundary 
point. An example of the computed tide compared to the actval 
tide is shown in Figure (7). 

Before actual tide data wer2 used as input to the 
meade, a test run was made to determine the length of running 
fame needed for the model to propagate the tidal wave through— 
out the entire bay. The M. (lunar-semi diurnal tide) sine 
wave recycles every 12 hours, 25 minutes. The Monterey M, 
consituent was used for the test runs. The tide travelled 
from the mouth of the bay to Moss Landing in less than a haif 
hour. This does not mean that all the perturbations caused 
by the initial wave placed on the grid have been dampened. 

In the original Hansen model the wave heights at all the gr:.d 
points were set equal to zero and then the initial wave was 
inserted at the open boundary points. It seemed that if the 
initial tide was at either high or low water, then large 
transients would be induced as the tide travelled across the 
grid. Several test runs were made to test if the entire grid 
should be initialized at the same tidal elevation rather than 
just at the mouth of the bay. It had been hoped that this 
would keep the initial shock of the induced wave to a minimum. 
Unfortunately, this technique caused a separate wave to be 
generated along the closed boundary. It was decided that for 
Monterey Bay all computer-runs would be made _so that the 


immtial tide would be at_mean sea level. Then, since the 


re —— _— 


tide moves quite rapidly across the bay, the transients 


would be kept to a minimum. 
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Figure 7. Computed vs. Actual Tide at Santa Cruz. 





fe OEEZOncal Eddy Diffusivity 
The parameter alpha (Equations 20-22) which is a 
function of the horizontal eddy diffusivity coefficient, 
mesh length and time step has been described as a horizontal. 
viscosity »Darameter (Laevastu and Stevens) and as a stability 
factor (Mungalls and Matthews Brcianseas 1OGo)eeeOrtazg (atter 
Lax and Richtmeyer, 1956) defines alpha as 


K, br 
qa-= 1-4 5 (25) 
DL . 





where, KL 1s the horizontal eddy diffusivity term. 
Oritz and Hansen claim that the numerical calculations 
converge analytically as alpha approaches 1, but are unstabie 
when alpha 1s equal to one. Alpha equal to one implies 
the eddy diffusivity coefficient Ky is equal to zero, i.e. no 
diffusion of momentum. Since the time step and mesh length 
are arbitrarily selected for each rroblem, alpha depends 
primarily on the variability of the horizontal eddy diffusivity 
term. A general solution does not exist for eddy diffusivity. 
Empirical solutions must be developed for specific areas of 
interest. Garcia (1971) claims that the value of horizontal 
eddy diffusivity for Monterey Bay ranges from 04. to aii 
dynes/cm*. The latter produces an alpha of less than 0.5 
which is unrealistic. The former value produces an alpha of 
approximately 0.94 which is far more ame eEOe, . Test runs 
were made with alpha equal to 0.95, 0.990 and 0.998 for the 
M, constituent and 0.900, 0.940 and 0.950 for the actual 


tide. The best results for the tidal heights appeared to be 
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obtained with alpha equal to 0.950. This corresponds to a 


value of K, equal to slightly less chan 107 dynes /cm* 


nh 


Current speeds were not satisfactory but if alpha was 
pemeecased, then the tidal heights were over-dampened. 


5. Bottom Stress 


— 


As previously mentioned the quadratic expression 
(Equation 12) used to represent bottom stress has been in 
existence for sometime; however, there seems to be some 
discussion as to what value should be assigned to the bottom 
stress coefficient (r). The bottom stress coefficient is 


defined as 


ater c (26) 


7 


where C (cm 2 /sec) is the DeChezy coefficient and is a 


BunetLlOn Of bottom roughness and bottom material. 
Experience in the coastal waters around the Netherlands 


(Dronkers, 1964) has shown that a reasonable range for r 


3 5 


nee oom 408-5 ole 10) with 2.7 x 10 > the most 


common value used. Although Hansen (1966) used the values 


3 3 3 


ei .8 x 10 pues aX 10 ~ and 8.6 x 10 ~ for r without comment, 


it appeared that 2.8 x ome was used when the model depth 


was equal to or less than 50 meters (25 fathoms). Mungalls, 


5 


et.al., used 4 x 10 ~ for Cook Inlet, Alaska which has 


Numerous mud flats and shoals. Pekeris and Accad (1969) 


Suggested that 2 x ho = should be used in shallow water. 


3 5 


Test runs were made with r equal to 3.0 x 10 a Omak Lom 


ie 3.2 xX Oe for Monterey Bay. It will be shown in the 
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following section that the errors in the current speeds 
calculated by the mathematical model appeared to be quite 
large and the adjusting of the bottom friction coefficient 
Gees not seem to change the results significantly. The 
menue of 8.6 x ite seems to be the most satisfactory value 
eO use. f 
6. Wind Stress 

ii ite le a2s"*knewn about bottom stress, less is 
known about. the interaction of the atmosphere and sea water. 
It is exceedingly difficult to make the micro-measurements 
needed to accurately determine the wind stress coefficient. 
meewind stress or drag coefficient (A) is defined similarly 
to the bottom stress coefficient. Munk and others used the 
herve of 2.6 x ee (dimensionless) for wind velocities from 
12 to 40 knots (Dronkers). The drag coefficient is a function 
of the elevations at which the velocities were measured. 
Based on Cardone's investigation (1970), @ drag coefficient 
oni xX Moe was selected for Monterey Bay. The bottom stress 
term originally included the density of sea water and the 
wind stress term included the density of air. Thus, when 
equation (1) was divided by the density of sea water, density 
was eliminated from the bottom stress term but not from the 


wind stress term. 
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TILT. COMPUTER MODEL RESULTS AND DATA ANALYSIS 


A. DATA ANALYSIS 

Tides, currents and winds were measured on 12, 13, 18 and 
19 August and 6-13 November 1970 for the purpose of evaluating 
the mathematical model. The tides were measured at the three 
stations mentioned in the previous section by mechanical gauges. 


Seeeencar time Crrors’ can occur on the tide records if the 





a 


Meeaanical gauges are nolL maintained properly. The Monterey 


~ 





and Santa Cruz tide records were timed by chronometers. The 
August tide records for Moss Landing were not accurately 
timed and contain errors. Thus, only the tidal amplitudes 
from the August records of Moss Landing were used for compar- 
ison. The tides were mixed (almost diurnal) during the 12- 
13 August time period, semi-diurnal]. during the 18-19 August 
time period and shifted from semi-cliurnal to mixed during tne 
6-13 November 1970 time period. 

Wind data obtained from land stations and the USNS 
De Steiguer are listed in Table A-2. The land stations were 
located at the shore line. The exact location of the USN 
De Steiguer at any given time can be determined by noting when 
the current meters were emplanted and the current drogues were 
launched. Since Monterey Bay is a small geographical area 
for large-scale climatological parameters and the wind records 
from the various stations around the bay appear to be in fairly 


good agreement, the complete wind records from Moss Landing 
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fieagure 3) were considered to be representative of the 
entire grid for each analysis period. 

Currents were measured in a Lagrangian sense by current 
drogues and, during portions of the survey, in an Eulerian 
Sense by current meters. The structure and application of 
the current drogues are described by Stoddard (1971). Drogves 
were not placed north of the Monterey Canyon as there were no 
facilities to track them. The tracks of the current drogues 
are shown in Figures (A-1ll) - (A-16) and the hourly speeds 
and directions are shown in Table 3. The location and depth 
of the current meters are shown in Figure 3. The relationship 
between the currents derived from the drogues and some of the 
current meters, winds and tides are show in Figures (Al7) - 
(A-43). Data from other current meters did not appear to be 
realistic and were not used. The currents are plotted on the 
tidal curves. The heads of the current and wind vectors are 
set in the direction in which the parameters are moving. 

Tidal currents can either be the reversing type where tre 
currents are zero at the anti-nodes and a maximum at the nodal 
points or the elliptical type where the current direction 
passes through all points of the compass and the speed is 
seldom zero. Reversing currents can be found along the coast 
or in rivers. Elliptical or rotary currents, as shown in 
Figure (8), occur in the open ocean or large embayments. The 
elliptical currents move in a clocx-wise direction in the 
Northern Hemisphere. Maximum speeds usually occur mid-way 
between turning points and minimum speeds at high and low 


waters. If surface currents are primarily induced by wind 
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Current meter data (Station 1 in Figure 3), plotted ina 
Lagrangian format (Format A-44) indicated that the currents 
at this location, with only slight deviations, flowed ina 
steady northwesterly direction. Thus, the NPGS data seems to 
andicate that the oceanic currents had more influence during 
these surveys on the circulation in Monterey than the tides 
and winds. 

Professor William Broenkow, Moss Landing Laboratories, 
has measured surface currents in the vicinity of Moss Landing 
with drift poles (personal communication). The measurements 
were made during April, May, November and December, 1970. ‘The 
memes were never tracked for longer than 3 hours and 15 
minutes. The maximum current speed occurred during flood tide 
and was perpendicular to the coast. Currents with speed up 
to 1.5 knots and moving parallel to the coast were measured 


meme OL the canyon. South of the canyon, the speeds never 


exeeded 0.7 knots. The maximum current speeds occurred 
Sweenet GUrang flood or ebb tide in some cases. In other 
cases, the maximum speed occurred at high water. There were 


mee, enough Observations at this point to accurrately correlate 
the currents near Moss Landing with the tides; however, it 
did appear that the current flowed southwest during ebb tide 
and northwest during flood tide. In most cases, the poles 
travelled along equal depth lines 

McKay (1970), towed a Geomagnetic Electrokinetograph (GEK) 
from a moving ship in Monterey Canyon during July and August 
1970 to measure the Currents. The GEK only performed satis- 


factorily in areas where the water depth was greater than 75 
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fathoms. McKay believed that the measured currents which 
ranged in speeds from 0.05 to 1.35 knots, were dependent on 
tides. Closer observation of the data indicated that the GEK 
was probably meaSuring extremities or eddies of the California 
or Beewelling current. In one case, which McKay stated that 
the currents were a result of flood tides, the currents flowed 
in a direction that would not be usually associated with 

flood tides. The average speed of the current was 0.5 knots. 
As with the other current measurements, most of the GEK 
measurements tended to follow the contour lines of the canyon. 

Defant (1961) indicates that if a strong current or up- 
welling exists in the adjacent oceén, standing vortices can 
occur in coastal bays such as Monterey Bay. The flow of the 
standing vortex is such that the currents on the seaward 
side of the bay follows the main current while on the landward 
Side the current iS opposite the méin current. Theoretically, 
standing vortices would have the seme water mass circulating 
within it and there would be no water transfer from the main 
Glierent. in actuality, this is not the case. Variations lin 
the oceanic current or upwelling will perturb the stationary 
vortex and renew the water circulating in the bay. 

Carter and Kazmierczak claims that a permanent closed 
circulation system which moves in a counter-clockwise, direc- 
tion extends from Ano Nuevo, north of Monterey Bay, to Moss 
Landing (Figure A-16). The average speed of the system is 
0.1 knots. Coastal engineers at Pacific Gas and Electric 
Company (P.G. & E.) believe that they have accumulated enough 


field data to confirm this phenomena (personal communication). 
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Over a three year period P.G. & E. has launched a number of 
drogues at the 6 fathom line north of Monterey Canyon. The 
drogues were launched during all three climatic seasons and 
at different phases of the tidal cycle. The drogues consis-~ 
tently follow the contour line towa:xds the northwest. The maxi- 
mum speed of the drogues never exceeded O62 KnOES and usually 
was quite less. 

ine field data collected by NPGS in the bay during August 
and November 1970 only implied, but did not absolutely prove, 
that the oceanic currents and upwelling were the primary 
forces of the circulation of Monterey Bay. It should be noted 
that the oceanic currents and upwelling were the primary 
AQriving forces of the circulation of Monterey Bay. The NPGE& 
data were collected during the transitory times between the 
upwelling and California currents and the California and 
Davidson currents respectively. To properly correlate the 
currents in the bay with the oceanic circulation, measuremer.ts 
were needed in the adjacent ocean. Twelve hydrocasts were 
made in the open ocean during 8-9 November and again during 
13-14 November 1970. Five current drogues were also launched 
on 8 November 1970. The results from this data are shown 
in Figures (A-45) and (A-46). The dynamic height calculations 
m@euethie first time period indicate that the oceanic current 
diverged, with speeds ranging from 0.1 to 0.6 knots, as it 
approached the coast. The drogues which were launched north 
of the canyon continued to move north fairly rapidly. The 
second set of dynamic height calculations showed a more 


complicated current pattern in the ocean. A large eddy was 
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present southwest of Monterey Bay. The current speeds for 

this set of data were Similar to that of the first set. The 
oceanic flow may be affected by Monterey Canyon. The current 
meee OL th] canyon moves north and vice-versa. The observa- 
tionsin the bay seems to substantiate the oceanic data. How- 
ever, there was very little erent data available for the 
north basin of Monterey Bay and more surveys need to be made 

in this area. Ideally, simultaneous current measurements 
should be made in the ocean and the bay on a seasonal basis 

to really understand the circulation patterns in the bay and 
the adjacent ocean. It is obvious that measurements in the 

bay without corresponding measurements in the ocean are really 
not significant. It should also be noted that the Hansen 
hydrodynamical model is not designed to handle density currents 
and their presence will not appear in the computer calculations 
Ber the bay. 

One point should be made about the tendency of the current 
measurements to follow lines of constant bathymety. This 
implies that potential vorticity is being conserved. Potential 
vorticity is comprised of several components. The planetary 
vorticity, or Coriolis parameter (f), is related to the 
rotation of the earth and varies as the flow of water moves 
from the poles to the equator and vice-versa. Ina small 
area such as Monterey Bay, planetary vorticity can be con- 
Sidered constant. Relative vorticity (%) is a function of 
the motion of a fluid relative to the earth and can be defined 


as 


Q> 
Q> 


u u 
C= oyt - (27) 
1 2 


Q> 
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The sum of the relative and planetary vortices is the 
absolute vorticity. The conservation of vorticy equation 


(Stommel, i966) is 


fetes 7 
=e = 0 (28) 


d_ 
dt 
where D is the thickness of the fluid. 


When equation (28) is integrated, it can be shown that 


potential vorticity is 


f+ 2 
D 





= conscant (29) 


If flow is along lines of constant bathymetry (D = con- 
stant), and f is assumed constant, then absolute vorticity, 


mets constant. 


Bee COMPUTER MODEL RESULTS 

The computer results were compared with the field data. 
The parameters compared were water elevation heights and 
phases and current speeds and directions. The results were 
fetced in two different formats by an automated plotter. 
One format shows the current circulation in the entire bay 
Memalyegiven time. Only every third row is plotted in this 
format and some of the calculated details, especially along 
the coast, are now shown. The other format shows the currents 
plotted on the tidal curve for individual grid points. The 
calculated water heights and phases compared quite favorably 
with the real tides at Monterey Harbor, Santa Cruz and Moss 
Landing as can be seen when Figures (A-47) through (A-52) are 


compared to Figures (A-3), (A-4), (A-7) and (A-8) and Table 
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A-1. The model indicated that the tidal range was less in 
the canyon than it was in shallow water as seen in Figures 
(A-53) and (A-54). Although this could be expected, it would 
have been more satisfying if there had been actual -tide data 
available in the canyon to indicate if this was correct. 

The calculated directions also seemed to be accurate, Some 
of the calculated circulation patterns in Monterey Bay are 
shown in Figures (A-55) through (A-53). As anticipated, the 
Meael produced two gyres ~ one in the northern basin and the 
other in the Seen basin - separated by the canyon. The 
gyres were more predominant during the ebb and flood periods: 
The northern gyre rotated in a counter-clockwise direction 
(Figure A-61) and the southern gyre rotated in a clockwise 
direction (Figure A-60) on the flocd tide and vice versa on 
the ebb tide. The maximum currents occurred at high water 
(Figure A-62) and low water (Figure A-59). At high water the 
currents diverged from the canyon into the basins and at low 
Watch tme CGUrErents Converged towarcs the canyon. The currents 
in the canyon are shown to be reversing currents. Although 
not shown in the figures, the model also depicted the currents 
along the coast-line as reverSing currents. 

With few exceptions, the calculated circulation pattern 
femadicimiliar to the circulation pattern depicted by the field 
measurements. This is another indication that the tidal 
current is dominated by the oceanic currents. 

The calculated current speeds appear to be an order of 
magnitude too large. Except for near Moss Landing, the maxi- 


mum current speeds measured in the bay was 0.7 knots and the 
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average speed was 0.3 knots. Since it has been concluded 
that the oceanic currents have more influence in Monterey 

Bay than the tides and winds, then it would be expected that 
the tidal currents would be in the order of 0.1 knots or less. 
There are two other reasons to expect that tidal currents 
would be low. First, the fact that the tidal phases are 
nearly equal throughout the bay implies that the horizontal 
current speeds would be low. Second, the magnitude of the 
current is a function of the water nass displaced horizontally 
during the rise and fall of the tides and the cross-sectional 
area through which it moves, (Simmons, 1966). Mathematically, 


this can be represented by 


(@.) th.) (£,) = AL (n,/*) (30) 


e 


where, u_is the average velocity at the entrance of the bay, 
h_ is the average depth across the entrance of 
Monterey Bay (approximately 612 feet), 
i lomenomrengtieon the Entrance of Monterey Bay 
(19.5 miles), 
Ay is the area of Monterey Bay (approximately 175.5 
square miles) 
(n,/T) is the average increase in water elevation in the 
bay during the time period T 
As an example, on 7 November 1970 there was an increase 
of two feet in water elevation from low to high water in 5 
hours and 30 minutes. Thenu. is approximately 0.05 knots. 


e 


For the same time period, the model calculated an average 
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speed 0.1 knots for the south end of the opening; 0.2 knots 
for the canyon; and close to 1.0 knots at the north end of 
the opening. Thus, it appears that the model does not 
correctly calculate the current speeds in Monterey Bay. 

The model computer program was checked thoroughly for 
errors. The continuity equation (Equation 11) and the con- 
servation of momentum equations (Equations 12 and 13) were 
checked by hand calculations. Decreasing the viscosity 
parameter (a) and increases the stress coefficient (r) dampens 
the current velocities; however, i* either term is varied 
TOO greatly, then the water elevations are over-dampened. 
This occurred in the test cases when a and r were set equal 
we 0.90 and 3.2 x roe respectively. These two parameters 
are useful only when the computer results need a fine 
adjustment. 

Durine most of the test runs, the computer program was 
only run out to 24 hours. It was thought possibly that even 
though all the tidal heights on the grid reached their 
correct height in a short time, the transients created by 
the initial waves moving across the grid would need more time 
to be completely damped. One pacer commencing 1900 PST, 

6 November 1970, was run out to 64 hours. The length of the 
computer run did not improve the calculations, which implies 
equilibrium conditions are obtained after about one hour. 

Several test runs were made so that the computer calcula-~ 
tions at every time step between 0 and 120 seconds and 7200 
to 7320 seconds could be observed. It became apparent that 


the velocities increased greatly as the waves moved from the 


Si 








a 


deep canyon areas to the much more shallow basins. Some 
mrregularities in the currents occurred along the very shallow 
(less than 5 fathoms) coastal boundaries. It was decided to 


smooth the bottom by the following technique: 


h ( = 0.5 One Zo. Ch + 


N,M) nM) eee ae 


(31) 


+ 


eee Ce) 


There was approximately a ten percent improvement in the 
current speeds, particularily Beet the canyon, when the 
bathymetry was smoothed once and approximately 15 percent 
improvement when the bathymetry was smoothed twice using 
eguation (31). It would appear that the model, employing tae 
grid spacing used here, has difficulty coping with the great 
variance in depth over as small an area as Monterey Bay. The 
large lateral and transverse slopes of the canyon (Figure 4 
amomuigure 5) are Significant bathymetric problems. It is 
believed that if a finer grid was used, the accuracy of the 
calculated current velocities would be improved. This point 
was proved in reverse when the grid was changed to 12 x 25 
with a mesh length of one nautical mile, the currents became 
larger. 

Assuming that the present grid 1s inadequate for Monterey 
Bay, some mention about the computer requirements for this 
model need to be made. Remembering that the grid size is 19 
x 46 and the time step is 6 seconds, the Control Data 


Computer (CDC) 6500 computer at FNWC utilizes 53000, computer 


8 


words, requiring 3.25 hours to compute a real time interval 
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Bre24 hours. £f a finer grid and, thus, a smaller time step 
was used, the computer core storage and running time would 
be greatly increased and the model would be impractical to 
use. 

At the suggestion of Dr. Tiavo Laevastu, FNWC research 
oceanographer, the bathymetry of all grid points which were 
less than 10 fathoms, were increased to 10 fathoms. This 
Suggestion seemed to be appropriate since the current speeds 
would then become less along the coast and possibly the lower 
speeds would be diffused into the rest of the bay. A test 
run was made for 24 hours. The calculated current speeds 
did improve at the head of the canyon and other areas where 
there was a large slope in the bathymetry, but in the flat 
basins the currents speeds became larger. 

A more serious problem than the large variation in the 
Monterey Bay bathymetry is the inability of the model to 
correctly calculate the currents along the open boundary of 
the grid. Equation (30) and analysed field data indicate 
that for Monterey Bay the average tidal currents across the 
open boundary should be less than C.1 knots. The current 
speeds calculated by the model were an order of magnitude 
larger. The difference between the calculated and actual 
current speeds seem to remain fairly constant throughout the 
Day. The inability to calculate the proper currents along 
the open boundary may be an inherent trait of the Hansen 
model. This author has found that the model produced 
unusually large currents at the entrances of the Gulf of 


Mexico and the South China Sea. Maury Pelto, National Marine 
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Fisheries Service, has observed the same phenomena when using 
the model for Bristol Bay, Alaska. The large currents at the 
open boundaries of these three models did not seem to be 
transmitted across the grid. However these areas were much 
Barger than Monterey BEY and the mesh lengths were of the 
Same order of magnitude as the enti:e Monterey grid. It may 
be that for larger geographical areas, the frictional terms 
in Equation (7) may dampen out the higher speeds as the tides 
move across the grid. 

There may be an additional problem with the Monterey Bay 
open boundary. The assumption that the tidal ranges and 
phases were: equal for all the grid points along the open 
boundary may be part of the reason that the current speeds 
were too large. As mentioned previously, since the Monterey 
and Santa Cruz tide records were in good agreement, the Santa 
Cruz tides were used along the open boundary. In most 
embayments, this would have posed no problems; however, since 
the canyon crosses the entrance of Monterey Bay, the tidal 
Benge im the canyon could be different than the range at the 
more shallow end points. If the tidal range was less in the 
canyon area, then the current speeds would have been smaller 
as the wave moved into the basins. It then becomes quite 
apparent that the tides and currents along the open boundary 
have to be more accurately defined. 

Since the calculated current velocities did not appear 
=EembDencOrrece, 1b did not seem apropos to run the model with 


wind fields. Only one twenty-four hour run was made with a 
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wind field included. The run began 1200 PST, 12 August 1970. 
The wind field had a velocity of 10 knots and a direction of 
225°. Altnough it is difficult to make a general statement 
about the influence of a wind field on Memeeney Bay based 

on one comouter run, in this case the tidal range seemed to 

be increased by approximately one foot and the current speeds 
mao tO 0.2 knots. It should be remembered that the actual 
tides and not the predicted tides are inserted at the open 
boundary and the wind effects are already included. By adding 
a wind field, the additional increases in water elevation 


and current speed result in slight over predictions. 
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IV. CONCLUSION 


As the field data analysis and the computer test runs 
progressed, it became quite obvious that the assumption that 
Monterey Bay would be an excellent site to evaluate the Hansen 
numerical model was invalid. First, field data collected by 
NPGS and other groups indicated that the principal driving 
force of the circulation in Monterey Bay was the adjacent 
oceanic current rather than the tid2s. At least one standirg 
vortex appears to exist in Monterey Bay. Except for the Moss 
Landing area, maximum current speeds in the southern basin 
me O./0 knots. Currents at the head of Monterey Canyon aré 
undoubtedly influenced by the tides where the currents reached 
Medsured speeds of 1.5 to 2.2 knots. The canyon may also 
Pmttuence the oceanic flow but there was not enough data to 
really confirm this hypothesis. Comprehensive field surveys 
need to be made on a seaSonal basis in the bay and adjacent 
ocean to PU vammderstand Fhe Ccireulation patterns of the 
area. 

If the density currents are the predominent currents in 
Monterey Bay, then the tidal-induced currents would be less 
Pidmeene waensity Currents and, except for the shallow coast- 
line, should be on the order of 0.1 knots. The computed 
velocities were as much as an order of magnitude too high. 

It is speculated that the steep slopes of the canyon and 


conditions imposed at the open boundary are the principle 
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mause of the large errors in the current speeds. A finer 

grid system may solve the problem but would not be practical 
on present-day computers. Smoothing the bathymetry and 
adjusting the horizontal eddy viscosity parameter and the 
meeconm Stress coefficient did not improve the current predic- 
tion results appreciably. Miereasing the shallow water depths 
to a minimum of ten fathoms did improve the results. The 
calculations of water elevation and phase and current direc-: 
tion which indicates that separate gyres exist in the northern 
and southern basens appears to be correct. 

This thesis points out the important fact that when the 
Hansen model or any other numerical model isS applied to 
Specific azeas cf interest, the calculations should be 
correlated with actual field data. Parameters such as the 
horizontal eddy diffusivity term and bottom stress coeffi- 
cient can be considered empirical in nature and vary from 
location to location. It is very important to properly define 
the tide and current input along the open boundary. If, 
aS in Monterey Bay, there is a large variance in the depths 
@eross the open boundary, then the mee should be accurately 
measured at all locations. Depending on the size and com- 
plexity of the area of interest, at least one and preferably 
several tide gauges should be used along the closed boundary. 
In the case of Monterey Bay, there was a noticeable lack of 
Mids edbame@qmtne Nortn basin. Thus, it was difficult to 
make definite conclusions about the circulation in this area. 


Finally, this report demonstrated the need to understand 


all the forces that govern the circulation in embayments. In 


57 





cases Similar to Monterey Bay, where either permanent or 
seasonal oceanic currents exists, one must look at the embay- 
Ment aS palit of a larger system. Ideally, computer models 
should include all the factors governing the circulation in 
an embayment but the amount of computer time reguired for 
these more sophisticated models makes this suggestion 
impractical on present-day computers. Thus, in applying anv 
model to an area, there should be some prior knowledge as to 


which model would be more significant for the area. 
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APPENDIX A 
Monterey Bay Tide, Current and Wind Data 
for 12-13, 18-19 August and 6-13 November 1970 
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Figure A-14, Tracks of Drogues 34- Z. 
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Figure A-15. Tracks of Drogues 35-38. 
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Figure A-55. Calculated Currents in Monterey Bay at 
Higher Low Water plus 0.5 hours. 
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Figure A-56. Calculated Currents in Monterey Bay at 
Lower High Water minus 2 hours. 
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Figure A-57. Calculated Currents in Monterey Bay at Lower 
High Water plus 1 hour. 
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Figure A-58. Calculated Currents in Monterey Bay at 
Lower Low Water minus 2.5 hours. 
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Figure A-59, Calculated Currents in Monterey Bay 
at Lower Low Water. 
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Figure A-60. Calculated Currents in Monterey Bay 
at Lower Low Water plus 3.5 hours. 


tele 





MONTEREY BAY 





oy Cll SeC 


— 


1.2@5 KNOTS 


Figure A-61. Calculated Currents in Monterey Bay 
at Higher High Water minus 2.5 hours. 
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Figure A-62. Calculated Currents in Monterey Bay 
at Higher High Water plus 0.5 hours. 
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Figure A-63. Calculated Currents in Monterey Bay 
at Lower-Low Water minus 2.95 hours. 
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Table A-I. Comparison of Tides at Monterey Bay Reference Stations 


|____ Monterey Santa Cruz Moss Landing 


Day/Time} High- |} Diff- |Day/Timel High- | Diff- pPay/Time | High- | Diff- 
(PST) erence (PST) | Low (PST) Low erence 




























erence 
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12/1930 12 8.55 a 
13/0300 I 2.8 i“ 
13/1000 13. 6.35 oe 
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18/1700 18 7.55 = 
18/2300 1S Soke _— 
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19/1200 19 2.65 ia 
19/1800 aces _— 
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Comparison of Tides at Monterey Bay Reference Stations (Cont'd) 































Monterey Santa Cruz Moss Landing 
Day/Time] High- | Diff- |Day/Time] High- | Diff- ! | Diff- 
(PST) | Low erence | (PST) | Low ; erence erence 
November | November November 
06/0503 6.6 
0.8 
06/1000 6/000" “Sad 
ks 
| 06/1600 067 1530 © 7 215 
4.55 
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Cm 
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| S50 
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4. 4.75 
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4, 4.45 
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| ae 353 
09/1900 09/1900 6.75 
334 
10/0100 110/0100 3.6 
4.7 
10/1400 10/1430 2.7 
) 3.95 
10/2000 1072000 6.65 
2./ 
11/0200 y0200 3.95 : 
4, 


124 





{ 


Day/Time 
GBS i) 


November 


11/0800 
11/1500 
11/2200 
12/0230 
12/1000 
12/1200 
12/2300 
13/0300 
13/1000 
13/1700 
14/0100 
14/0430 
14/1100 
14/1830 


Monterev 


14/1700 


Santa Cruz 





| Diff-  |Day/Time} High- 
November 
11/0800 =8.10 
fe) 
S00 sk. 35 
4.9 
Nee ZO .00 
cua 
12/0200 io 
4.8 
12/0830 ae 
19 
12/1600 Oo 
oo 
12/2200 
233 
13/0300 
a5 
13/0900 
Se 
13/1630 
as 
Wy 25919) 
{|e 
1/14/0300 
525 
14/1000 
Times! 


£25 


Moss Landing 


Day/Time} High- 
(PST) | Low 


November 


11/0800 
11/1500 


11/2100 
1 12/0200 


12/0800 
Zo ae 
12/2230 
13/0300 
13/0900 
13/1630 
13/2300 
14/0330 


11471000 


14/1700 


ome 


—!) 


Oo FSP DN 


5 
ao 
BD 


3 


5 


Deas 


cS 
co 


70 


Comparison of Tides at Monterey Bay Reference Stations (Cont'd) 
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Time 


beogue 1: 


T2300 
1900 
2000 
2100 
2200 
2300 

130000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 


Drogue #2: 


122000 
2100 
2200 
2300 

130000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 


Drogue #3: No track (never gained contact) 


Drogue #4: 


T0200 


130700 


Aug. 


Aug. 


Aug. 


Table A-III 
DROGUE COURSE/SPEED DATA 


G7 


143 


1356 
160 
182 
ea 


sD 


Approx. Course Approx. Speed (kts) 


OZi5 
Oa) 
d0> 


0520 
Oe 210) 
O20 
O20 
0.30 
0220 
OG 
ORIZ 5 


O20 
0.34 
0.05 
OF EO 
O20 


0.15 
OnilS 
0.10 
0.10 
Vet > 
0.10 





0800 
0900 
1000 


Drogue #5: 


130200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 


Drogue #6: 


LSZH 
2200 
2300 

190000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 
oe 
1200 
1300 
1400 
1500 


Drogue #7: 


182100 
2200 
2300 

190000 
0100 
0200 
0300 
0400 
0500 

190600 
0700 
0800 
0900 
1000 


Aus. 


Aug. 


Aug. 


026 
028 
030 


2 
169 


168 
es 


094 
097 


eZ 
090 
124 
238 
Zak 
215 
129 
132 
043 
O97 
065 
056 
069 
130 
134 
144 
143 
i 
3 


Oy 


334 
Sa 
342 


Oi 
034 
034 
034 


Oe 


107 
163 


136 


Ono 


0.15 


0.40 
0.40 
0.34 
rea) 
OieES 
Om 
OerZ0 
OF 20 


OreZ > 
O25 
Or 0 
ORT 5 
Ora 
O25 
0.075 
Om 2S 
0.10 
O20 75 
O2025 
0. 30 
02220 
O30 
O25 
O23) > 
Om Zo 
Om225 
Omz7 3 


0225 
Oa 
0.20 
O20 
0.30 
D2 7S 
0.40 
O2272 
Omae 
0.425 
0.50 
0.425 
Oe 
OaZo 





1100 200 0225 
1200 096 02225 
1300 081 OAS 
1400 084 Oe275 
1500 143 O20 
Drogues #8: 
182200 Aug. 034 OZ 0 
2300 ae 0325 
190000 328 O42 5 
0100 328 Onis 
0200 003 0.45 
0300 030 0:20 
0400 071 0.20 
0500 089 On275 
0600 080 Or 
0700 087 OraZ > 
0800 090 025 
Meogue #9: No track (lost contact after 15 min.) 
Drogues #10: 
182300 Aug. 214 OS ALIAS: 
190000 252 0.40 
0100 285 0.175 
0200 356 O30 
0300 003 0.45 
0400 Ey OS 
0500 065 O27 5 
0600 093 On 17> 
0700 087 ORS 
0800 146 0.20 
0900 165 O25 
1000 180 ORS 75 
1100 207 0.40 
1200 235 O32 5 
1300 237 0.225 
1400 263 0.05 


Drogues #1ll: Very short track. Moved due south 250 yds., 
2500-2400. Sshiftvceaseco approximate head 250 
for 30 min. covering 200 yds. Reversed track 
Eomapout O70 Lom Zens. covering 400 yds. 


Lost contact at 190230 (shadow zone?) 


Direauer fila: 
190000 Aug. 230 Oe2 75 
0100 229 02525 
0200 261 0.425 
0300 293 O2aS 
0400 026 0.40 
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Drogue #13: 


Drogue #28: 


0500 
0600 
0700 
0800 


190000 


0100 
0200 
0300 
0400 
0500 
0600 
0700 


061300 


1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
23010 


070000 


0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 


180000 


Aug. 


Nov. 


056 
012 
064 
ly 


DOIo 
257 
305 
DZ 
028 


016 
068 


Saw) 
Ban 


311 
Biz” 
3277 


Zouk 
254 
244 


248 
255 
205 
280 
273 
298 
315 
327 
330 
352 
348 
eo 
342 
ag 0 
S2y) 
342 
342 
336 
2a 
250 
248 
247 
247 
247 
276 
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O2425 
Ue 
OR 

O34175 





Drogue #29: 


Drogue #30: 


Drogue #31: 


0100 
0200 
0300 
0400 
0500 
0600 


061200 
061300 
1400 
Pol 
1600 
1700 


071800 


1600 


071900 
2000 
2100 
2200 
2300 

080000 
0100 
0200 
0300 
0400 


Nov. 


Nov. 


Nov. 


Zo.) 
215) 
292 
314 
316 
304 


048 
048 
OSe: 
040 
056 
082 


268 
265 
244 
33 
207 
241 
249 
248 


23 
263 
one 
302 
345 
obe 
Baie 
002 
002 
348 
310 
276 
264 
244 


29:3 
279 
278 
Ze 
221. 


290 
212 
SUS 
294 
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OOOO OC O 
OJ 
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OOOO CO © 
J 
Sa) 


ooooooo°oo°oo°o°o°oeooo0o0oo00 80 O 
NO 
LS 





0500 SIL, OZ 2 


0600 339 0.65 
0700 002 0) AS 
0800 010 0.34 
0900 Or. 0.475 
1000 046 0.40 
1100 ogl O25 
081200 009 Orly > 
1300 a3 Oe 
1400 335 0.10 
1500 250 Oe Aro 
1600 So8 Oe 25 


Drogue #32: Very short track. Moved for 1-1/2 hours on 
approximate heading 170 covering 550 yds. 
contact (shadow zone?) 


Drogue #33: 


072000 Nov. 095 ORaEyS 
2100 7 O20 
2200 ee 0215 
2300 229 O2325 

080000 234 Orz7 
0100 Zar, 0.34 
0200 Ze 0530 
0300 284 02475 
0400 S25 DE 2S 
0500 344 D225 
0600 359 OE 
0700 012 OFS 15 
0800 030 On225 
0900 048 0.40 
1000 Oral O25 
1eOID 043 Oe 2e7S 
1200 Oa oO. Ie 
1300 041 O25 
1400 054 OnZ 75> 
1500 267 Orvis 

Drogue #34: 

101000 Nov. 246 0.34 
1100 23.3 0.45 
1200 238 Oe si 
1300 226 0.40 

Dregue: x 

100800 Nov. 240 0.25 
0900 234 0.34 
1000 Zion 152 
1100 Zoe Ore 7 
1200 Delp O45 
1300 244 0.34 
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Peogue 
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oye" e 


Drogue #35: 


Drogue #36: 


Drogue #37: 


Drogue #38: 


100900 
1000 
1100 
1200 
1600 
1700 
1800 
1260 
2000 


101100 
1200 
1600 
1700 
1800 
010 
2000 


132000 
0300 


0-300 
0400 
0500 
0600 
0700 
0800 


130500 
0600 
0700 


Nov. 


Nov. 


Nov. 


Nov. 


Nov. 


Nov. 


190 


203 


180 
164 
SZ 
io 
201 
223 
ZA. 


210 
ZS 
Z05 
Zag 
241 
210 


GZ 
180 
205 
214 
Ze 
214 


219 


228 
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0 
0 
0 


Litieeisistipted. Track 


O97. oo Oo © © © 2 o@ 
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oueD 
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nei 
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